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ABSTRACT 

We present ground-based optical and Spitzer Space Telescope infrared imaging 
observations of the ecliptic (Jupiter-family) comet 21P/Giacobini-Zinner, the 
parent body of the Draconid meteor stream, during its 2005 apparition. Onset of 
nucleus activity occurred at a pre-perihelion heliocentric distance, ~ 3.80 AU, 
while post-perihelion 21P was dusty (peak Afp = 131 cm -1 ) and active out to 
heliocentric distances >3.3 AU following a logarithmic slope with of —2.04. 
Coma colors, V - R = 0.524 ± 0.003, R - I = 0.487 ± 0.004 are redder than 
solar, yet comparable to colors derived for other Jupiter-family comets. A nucleus 
radius of 1.82±0.05 km is derived from photometry at quiescence. Spitzer images 
post-perihelion exhibit an extensive coma with a prominent dust tail, where 
excess emission (over the dust continuum) in the 4.5 pm IRAC image arises 
from volatile gaseous CO and/or C0 2 . No dust trail was detected (3cr surface 
brightness upper-limit of 0.3 MJy sr _1 pixel -1 ) along the projected velocity vector 
of comet 21P in the MIPS 24 pm image suggesting that the number density of 
trail particles is < 7 x 10~ n m -3 . The bolometric albedo of 21P derived from the 
contemporaneous optical and Spitzer observations is A{9 = 22°) = 0.11, slightly 
lower than values derived for other comets at the same phase angle. 
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Subject headings: Comets: individual (21P/Giacobini-Zinner): infrared: solar 
system: meteors, meteorides 



INTRODUCTION 



Comet nuclei formed beyond th e protoplanetary disk frost line (heliocentric distances, 
r h >5 AU; iLunine fc Gautierl 12004 ) . among the giant planets and were scattered into the 
Kuiper Belt and beyond into the Oort Cloud (OC). Since their formation, the interiors and 
surfaces of most comets have re mained at temperatures below 140 K while in "cold storage" 
in the Kuiper Belt or the OC (jMeech fc Svorehl 120041 ) . Moreover, most nucleus surfaces 
have remained below 400 K even during perihelion passage. At such low temperatures, dust 
mineralogy remains stable and each comet nucleus retains a record of the minerals, ices, and 
volatiles extant in the comet agglomeration zones in the early solar system. Furthermore, 
comet nuclei may retain their primordial compositional inh omogeneities so that dif ferent to- 
pographic regions would have different compositions (e.g.. lDello Russo et al.ll2007l ). leading 
to variations in volatile production rates as a function of nucleus rotation. Nucleus het- 



erogeneities are app arent in the 



distinct topography fjBelton et al. 



( IFeaga et al.ll2007l : lA'Hearn et al 



y-by imagery of ecliptic comet 9P/Tempel 1. Regions of 
2007), and spatially distinct sites of water and CO2 release 



subsurface composition (IHarker et al.l 120071 : iKadono et al.ll2007l ). 



200 5]) are evident, as well as h eterogeneities in surface and 



There are two general dynamical families of comets, classified by derived orbital ele- 
ments (i.e., Tisserand parameter, Tj) based on current observations. Nearly-isotropic comets 
(NICs; Tj < 2) have orbits that are approximately uniformly distributed in inclination, and 
are derived from the OC. Ecliptic comets (EC) have orbits that are confined to inclinations 
near the ecliptic plane, 2 < Tj < 3, and originate in the Kuiper Belt. Resulting from frequent 
perihelion passage over the 4.0 - 4.5 Gyr period since their formation, EC comets have be- 
come noticeably less active than OC comets, characterized by lower gas and dust production 
rates. Multi-epoch spectral energy distributions (SEDs) of ECs, from which comae dust 
properties can be constrained, are needed to assess the possible interrelationships between 
their reduced activity levels and dust properties. 

The study of the physical properties of cometary nuclei and comae both are equally 
important to our understanding of the outer solar system environment during the era of 
icy planetesimal formation, and complement efforts to discern conditions extant in early 
protoplanetary disks during the epoch of planetesimal formation. We present new optical and 
contemporaneous Spitzer observations of comet 21P/Giacobini-Zinner obtain ed during its 



2005 apparition obtained as part of a larger survey of both ECs and NICs (e.g.. iKelley et al. 
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20061 ). Ground-based optical observations enabled us to obtain precise optical photometry 



to asses variations in dust productivity with heliocentric distance and to study the the 
comet's near-nucleus structures, including its jets and coma. Spitzer images at mid-infrared 
wavelengths enable the study of the spatial distributions of volatiles and dust in the coma, 
as well as facilitating investigation of comet trail properties. We describe our observations 
in £j2j and discuss analysis of our optical and infrared (IR) observations in §3], while our 
conclusions are summarized in © 



2. OBSERVATIONS AND REDUCTION 

Comet 21P/Giacobini-Zinner was discovered by the French astronomer Michel Giacobini 
in 1900 and rediscovered two apparitions later, in 1913, by the German astronomer Ernst 
Zinner. 2 IP, the parent body of the Draconids (also know as the Giacobinids) meteor shower 



( lBeechlll986l ). is a Jupiter-family (i.e., EC with a short-period) comet with an orbital period 
of 6.61 years and an aphelion distance just exterior to Jupiter's orbit. In 1933 a spectacular 
Giacobini meteor storm was visible across Europe. The 1946 apparition of 21P was especially 
noteworthy as the comet passed only 0.26 AU from Earth in late September with a visible 
magnitude ~ 7. In early October, an unexpected outburst of activity caused the comet 
to brighten to sixth magnitude. Every alternate apparition of comet 2 IP presents favorable 
observing geometries for observers on Earth, with integrated coma visual magnitudes as high 
as 7. 2 IP was also the first comet visited by a spacecraft when the Int ernational Cometar y 



Explorer (ICE) flew past at distance of 7,800 km on 1985 September 11 (IBrandt et al.lll9881 ). 



2.1. Ground-based Optical 

Optical images of comet 21P/Giacobini-Zinner were obtained on the University of 
Hawai'i (UH) 2.2 m telescope on Mauna Kea, using a Tektronix 2048 x 2048 CCD camera 
at the f/10 focus of the telescope (image scale of 0.219" pixel -1 ) through a Kron-Cousins 
filters set (V: A = 5450A, AA = 836A; R: A = 6460A, AA = 1245A; /: A = 8260A, 
AA = 1888A), during multiple observing runs from 2004 June, 2005 October and December, 
and 2006 March. The camera read-noise was 6.0 e~ with a gain of 1.74 e~ per ADU. Non- 
sidereal guiding at the cometary rates of motion, derived from the JPL /Horizons ephemeris, 
was used during all image observations. Specific observational details, including orbital ge- 
ometry for 2 IP, filters, number of exposures, individual exposure times, and sky conditions 
are summarized in Table [TJ 
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The optical CCD frames were reduced with standard IRAF0 routines. All the images 
were reduced using flat field images taken during the evening and morning twilight sky and 
cleaned o f bad pixels and cosmic rays. The frames were calibrated with the standard stars of 



Landoltl (119921 ). which were observed on each photometric night. Observations of typically 



20 standard stars were obtained over a range of air-masses, and with a wide dispersion of 
color to fit for both extinction and color terms. Apparent magnitudes of twenty to thirty 
stars field stars of equal or greater brightness to the comet on each frame were measured in 
order to do relative photometry. After correcting the measured magnitudes for extinction, 
we used the deviations of the field star magnitudes in each frame from their nightly average 
values to correct for frame-to-frame extinction in the comet's measured signal. Figure [1] 
shows selected individual optical images for each ground-based observing run. 



2.2. Spitzer Imaging 



We observed comet 2 IP with the Infrared Array Ca mera (IRAC 



the Multiband I maging Photometer for Spitzer (MIPS; iRieke et al. 



azio et al.ll2004l ) and 



20041 ) on the Spitzer 



Space Telescope (IWerner et al.ll2004l ). post-perihelion as summarized in Table HJ The MIPS 
astronomical observation request (AOR) used the 24 /im array in mapping mode to assess the 
existence of the comet's dust trail (AOR key 0015734784, 22' x 15' map). A duplicate MIPS 
observation was taken 25 hr later to observe the sky background after the comet had moved 
out of the frame (AOR key 0015734528). The IRAC AOR utilized the 4.5 and 8.0 fim arrays 
to obtain complementary deep photometry of the coma (AOR key 0013820160, 6.7' x 5.5' 
map). 

The IRAC and MIPS observations were calibrated with Spitzer pipelines S13.2.0 and 
S13.0.1, respectively. The images were mosa icked in the rest frame of the comet with 
the MOPEX software jMakovoz & KahiJ boOa ) at the native IRAC and MIPS pixel scales 
(1.22" pixel" 1 for IRAC, 2.5" pixel" 1 for MIPS 24 /im). The MIPS background observation 
was mosaicked with the same size, scale, and orientation as the primary observation, then 
subtracted from the primary. Before mosaicking, all images were masked to remove cosmic 
rays and bad pixels. No additional corrections were necessary. 

Our IRAC image of comet 21P was obtained at a post- perihelion dis tance of = 
2.40 AU. Analysis of the secular optical light curves of 21P by iFerrinl (120051 ) suggests that 



1 IRAF is distributed by the National Optical Astronomy Observatory, which is operated by the Asso- 
ciation of Universities for Research in Astronomy (AURA) under cooperative agreement with the National 
Science Foundation. 
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the nucleus turn-on time (marking the onset of steady coma activity) occurs pre-perihelion 
near 3.7 AU, whereas onset of nucleus turn-off may not occur until 21P reaches 5.4 AU. Indeed 
Tancredi et aP fepQoh argue that 21P is active at 4.5 AU. The IRAC 4.5 pm image, Fig. [2J 
shows an extended coma surrounding the nucleus. Discussion of this extended emission is 
presented in §3.71 The Spitzer MIPS 24 /im image is presented in Fig. [3j 



3. DISCUSSION 
3.1. Optical Photometry 

The photometry derived from our optical images enable investigation o f comet ary ac- 



tivity with heliocentric distance, estimation of dust production rates (Afp; lA'Hearn et al. 



19841 ). determination of the coma color, and construction of rotational light curves. 



The stellar and comet fluxes were extracted from the optical CCD frames using the 
IRAF photometry routine "PHOT" with a circular aperture. This routine automatically 
finds the centroid of the image within the user-specified photometry aperture, with the sky 
background determined in an annulus lying immediately outside the photometry aperture 
(for stellar images), or selected from an average of many sky positions outside the extent 
of any coma. Photometric apertures between 1"- 5" were used for the comet photometry 
with the average sky background determined using sky annuli with an inner radii of 8.0" 
- 10", with widths of 4.0" - 5.0". Averaging many sky background positions adjacent to 
the comet also enabled rejection of any bad pixels or field stars found in individual sky 
annuli circumscribing the aperture used to determine the comet surface brightness. The 
smallest possible aperture that minimized contamination from extended surface brightness 
of the coma while including most of the flux from the nucleus (with respect to the full-width 
half-maximum seeing determined from stellar point sources on the frame) was 3". A 4" 
aperture was also used for the field stars used for differential photometry. After correcting 
the measured magnitudes for extinction, we used the deviations of the field stars from their 
nightly average magnitudes in each frame to correct for frame-to-frame extinction variations 
in the comet's measured signal. 

For the non-photometric data (Table [1]) from the nights of 2004 June 21 and 22, and 
2006 March 05, the comet fields were re- imaged under photometric conditions on 2005 May 
18 and 2006 June 03 for "boot-strap" calibration. Multiple ( > 20) field stars were measured 
on both the calibration images and the non-photometric images and differential photometry 
was used to compute the photometrically-calibrated brightness of the comet . This technique 
works well for up to ~ 0.5 mag of extinction. 
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3.2. Heliocentric Activity 

The most direct indicator of nucleus activity is the appearance of visible coma around a 
comet nucleus. Our first observations of 21P (Fig.QJi) were obtained on 2004 June (—375 days 
before perihelion) when the comet was at a large heliocentric distance, = 3.80 AU. No 
coma was evident on individual images, nor on the composite image created by mosaicking 
the 2 to 3 images obtained on each night to increase the signal-to-noise ratio of any diffuse 
coma emission. Absence of a det ectable c o ma — 375 days before perihelion is consistent with 



1991 apparition observations of iMuellerl (119921 ) who found no coma at —367 days before 



perihelion, yet a clearly detectable coma at —334 days. Our observations indicate that the 
heliocentric turn-on point for 2 IP nucleus activity this apparition was similar to the behavior 
observed in the 1991. 

The average i?-band magnitude measured for 2 IP on a given night are summarized 
in Table El column [3], as is the reduced magnitude, column [6]. The reduced magnitude, 
m(l, 1,0), i.e., the observed magnitude normalized to unit heliocentric distance, r h , and 
geocentric distance, A, for zero phase angle, a, was derived from the relationship for active 
comets 



m(l, 1, 0) = m R - 2.5 n log(r h ) - 2.5 k log(A) - a(3, (1) 

where tur is the i?-band photometric magnitude measured in a 3" circular aperture, (3 the 
linear phase coefficient = 0.035 mag deg -1 , and n and k are constants, taken to be unity 
and 2 respectively. Nominally, values for these two latter constants (n being dependent on 
dust grain production, coma shape and size) are derived from modeling a given comet at a 
variety of heliocentric distances if there are sufficient observational data. Equation [1] also is 
appropriate for expressing the brightness of an active comet wi th a spherical coma aro und 



the nucleus (measured in a circ ular aperture , as described by iMeech Svorenl 120041 ). as 



ap 
in 



opposed to the relation used by iFerrinl (120051 ) which is the proscription for the brightness 



of bare nuclei (non-active stage). Computation of the reduced magnitude is a standard 
method that enables of photometric data of an individual comet obtained at a variety of 
heliocentric distances (different activity level) to be intercompared. For 21P, we obtained 
observations of the comet when there was no visible coma, when the comet exhibited a small 
coma around the nucleus, and when the comet was highly active displaying a significant, 
spatially extended coma with an elongated tail. To properly intercompare variations in the 
integrated aperture surface brightness, the i?-band photometric magnitudes were determined 
by measuring the comet's brightness in a small 3" circular aperture centered on the nucleus. 
Use of this small nucleo-centric aperture minimized the coma contribution to the measured 
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flux density, and allowed estimation of the nucleus brightness consistent with the uniform 
azimuthally-averaged surface brightness profile. 

Our observations of 2 IP span a range of heliocentric distance, both pre- and post- 
perihelion (Tabled]), thus we can examine trends in m(l, 1, 0) with r^. Figure H] shows the 
variation in the reduced magnitude as a function of derived from .R-band photometry, 
where the dotted horizontal lines represent the likely brightness range for a bare nucleus, 
calculated from our photometry data when no coma was seen. Evidently, cometary activity 
commenced near ~ 3.80 AU this apparition of 21P, and the majority of our optical 
observations and all of our Spitzer images were obtained while the comet was in an active 
state. 



3.3. Coma Colors 



Average coma color differences (V — R, R — I) are also listed in Table [2] for dates 
when mult i- wavelength observations were conducted. The average colors of comet 2 IP, 
V — R — 0.524 ± 0.003, R - I = 0.487 ± 0.004 are redder than the respective solar colors, 
0.36 and 0.28, in the Kron s-Cousi n s syst em, transforming the Johnson solar colors ((Allen 
19731 ) using the relation of iFernie Band emission from gaseous molecular species, 

such as C 2 (A « 5520A), NH 2 (A » 6335A), and CN (A « 9180 A), or atomic O I (A = 6300 A 
X D) potentially could contribute significantly to the flux den sity of comet comae observed 
within our VRI broadband filters. However, optical spectra of iFink fc Hickd (119961 ) indicates 
that band and line emission from these species is relatively weak in the coma of 2 IP. 

2 IP also lies along th e (R — I) versus (V — R) color-color trend line of Jupiter-family 
comet nuclei attributed by ISnodgrass et al.l (120051 ) to albedo variation (in § 13.51 we derive an 
dust albedo estimate for 21P), although we cannot quote a precise value for the nucleus color 
of 21P due to coma contamina tion. Broad-band colors of cometary nuclei are related to sur- 
face properties of the nucleus (jDavidsson & Skorovl 120021 ) and provide a metric to compare 
cometary nuclei to other small solar system bodies and to assess which parent body popu- 
lations give rise to Jupiter-family comets. The bluer color of 21P and other Jupiter-family 
comets, as compared to the me an color of Kuiper-Belt Object/Trans- Neptunian Object pop- 
ulations, V — R = 0.60 ± 0.07 (ISnodgrass et al.ll2005l ). suggests that cometary activity and 
surface processing in the inner solar system have altere d the primitiye surfaces of these n uclei 



agglomerated from materials in the Kuiper Belt (e.g., I Jewittl 120021 ; IWooden et al.l 120041 ) 
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3.4. Dust Production Rates 



The quantity Afp of a comet, a proxy describing the dust production rate (lA'Hearn et al. 



19841 ). is often used to describe nucleus activity vigor, and to distinguish between dusty and 



inactive comet taxonomies. Afp is given by 

Afp = (2.467 x 10 19 ) (2) 

a F Q 

where A is the wavelength-dependent albedo of dust particles, / is the filling factor of grains 
(total cross-section), p is the linear radius of the aperture at the comet, r h and A are 
heliocentric and geocentric distance (AU), a is the angular diameter of the field of view (in 
arcseconds), F f, s is the observed cometary flux and F & is the solar flux. The Afp parameter 
is typically expressed in centimeters. By providing a metric to quantify a comet's level of dust 
output (large Afp values indicate higher activity), one can ascertain the dust contribution to 
a set of photometric measurements and assess the heliocentric dependence of dust production 
rates in a given comet. 

i?-band values of Afp for 21P derived from our photometry are summarized in column 
[7] of Table [2] and are plotted as a function of in Fig. [5j Peak dust production rates of 
131 cm -1 occurred just post-perihelion. We fit the post-perih elion A fp values w ith a function 



of the form Gr\ using the least absolute deviation method (IPress et al.lll992l ). Our best-fit 
values are C = 414.4 cm, k = —2.04. The pre-perihelion observations do not span a large 
enough range of to derive a slope. Our post-perihelion Afp slope (1.76 < ry,,(AU) 2.92) 



is comparable to the pre-perihelion value of —2.08 ± 0.15 derived by lA'Hearn et al.l (1 995) 



based on observations obtained when 21P was at a ~ 1.34 AU. However, lA'Hearn et al. 



( 119951 ) find a significantly larger value (k = +0.38 ±0.16) when the 21P is interior to 1.86 AU 



post-perihelion suggesting that the comet's activity peaks just after perihelion passage. 



3.5. IR Photometry, Optical to IR SED, and the Dust Albedo 



We measured the flux of comet 21P in the MI PS and IRAC im ages using the aper.pro 
procedure of the IDL Astronomy User's Library (jLandsmanl Il993l ). The 2005 December 
IRAC, MIPS, and optical photometry are all measured at different heliocentric distances 
(rh = 2.1 — 2.4 AU), observing geometries (A = 1.7 — 1.9 AU), and spatial resolutions. 
However, differences in observational circumstances are small such that simple scaling laws 
can be applied to the photometry to create a composite 0.5-24 /im SED of comet 21P. For 
reasons discussed below, we elect to correct the IR photometry (§ 12. 2p to match the observing 
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conditions of the optical measurements (§ 12. ip . 

First, we correct the IR photometry for the limited resolution of the Spitzer telescope, 
as compared to the optical data. The IRAC and MIPS images are diffraction limited and 
stellar sources have full-width half maximum « 1.4", 2.0", and 6" at 4.5, 8.0, and 24 pm 



Fazio et al. 


2004; 


Rieke et al. 


2004) 



to compare to the 3" i?-band photometry in Table [2] (column [3]), and the IRAC PSF 
is marginally appropriate. To derive the IR photometry in an aperture size of 3741 km 
(corresponding to the optical photometry, 3" at 1.72 AU), we measure the comet's IR coma 
with increasingly larger apertures and use a curve-of-growth analysis to estimate the fluxes 
at apertures comparable to or smaller than the resolution of the Spitzer instruments. We 
examined the logarithmic profile of the optical data to verify that the coma slope at moderate 
sized apertures (p ~ 10") extends into the inner-coma (p ~ 3"). The photometry at 10" and 
the coma logarithmic slopes are presented in Table [3j 

Next we correct the IRAC and MIPS photometry for the differences in observer-comet 
distance. The correction is simply (A opt /A/^) 2 : 0.93 for IRAC, and 1.18 for MIPS. 

Finally, we scale the IRAC and MIPS photometry to match the heliocentric distance 
(77, = 2.32 AU) of the optical observations. There are two effects that depend on r^: the 
temperature of the dust, and the dust production of the comet. To correct for changes 
in the dust production, we scale the photometry using the Afp heliocentric distance trend 
derived for 21P, r^ 2M ( §3.41) . The MIPS photometry is scaled by a factor of 0.86 and 
the IRAC photometry scaled by 1.08. To correct for the dust temperature, we scale the 
photometry by B\{T opt ) / ' B\(Tir) , where B\ is the Planck function, and T opt and Tjr are 
the dust temperatures at the optical and IR heliocentric distances, which we set to the 
blackbody temperature, T = 278 r^ ' 5 K. The scale factors are: 1.37 for A = 4.5 pm, 1.20 
for A = 8.0 pm, and 0.88 for A = 23.7 pm for the IRAC and MIPS data respectively. 

Our measured fluxes, azimuthally averaged surface brightness profile logarithmic slopes, 
and the corrected photometry is presented in Table [31 The total SED of the comet is 
presented in Fig. El The optical magnitudes of Table [2] were converted to flux units using the 
zero magnitude flux densities of 3600 Jy (V-band), 3090 Jy (_R-band), and 2550 Jy (/-band). 
We fit the SED with a scaled, reddened solar spectrum to represent the scattered light at the 
optical wavelengths, plus a scaled Planck function, to represent the thermal emission from the 
IR wavelengths. We exclude the gas contaminated 4.5 pm data point ( §3. 7[) from the fit. Solar 
spectrum reddening ( the spectral slope) is typically expressed as percent per 0.1 pm (e.g., 



Jewitt fc Meechlll986l ). A negative spectral slope represents a bluing of the spectrum. Our 



observations constrain the scattered light r e ddeni ng between the V- and -R-bands, and the 



R- and /-bands. Following I Jewitt fc Meechl (119861 ). we express the reddening with the linear 
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slope S"(Ai,A2) expressed as % per 0.1 /mi. At optical wavelengths, our computed slopes 
for comet 21P (see Tabled are S'(V, R) = (15 ± 1)% per 0.1 //m, and S'{R, I) = (11 ± 1)% 
per 0.1 /an. The former slope of the reflectivity gradient is i n good agreemen t with the 
average value of S'(B,R) = (13 ± 5)% per 0.1 /mi derived by IJewitt &: Meechl (119861 ) for 
comets, while the latter S'(R, I) slo pe value agrees with th e general trend of bluer sca t tering 
at longer wavelengths (see Fig. 4 of IJewitt fc Meechl 1 1 9861 ) . Indeed, IJewitt &: Meechl (119861 ) 
observed 21P in the near-IR and found S'(J,H) = 3.5% per 0.1 /im and S'(H,K) = 0.3% 
per 0.1 /im. The fit has no degrees of freedom, so we do not present errors or a \ 2 analysis. 
The best-fit dust temperature is 206 K (~ 15% warmer than a blac kbody). We derive a 
bolometric albedo for the dust of comet 21P using Eqs. 6 and 7 from lGehrz fc Neyl (119921 ). 
The bolometric albedo of 21P is A(8 = 22°) = 0.11, sligh tly lower than values d erived for 
other comets at the same phase angle (A(6 = 22°) m 0.16, iKolokolova et all 120051 ). but not 
unusual. 



3.6. Nucleus Size & Light Curves 

The size of the equivalent spherical nucleus, Rn, can be estimated from the apparent 
mean magnitude of comet 21P measured at heliocentric distances when the comet is relatively 
inactive. Pre-perihelion, rriR = 21.89 ± 0.04 (Table [2]) and Rjy can be derived from the 
standard relation 



A R R 2 N = (2.238 x i 22 ) r 2 A 2 lO a4 [ m 0- m « +/3Q ] 



(3) 



where Ar is the geometric albedo, (AU) the heliocentric distance, A (AU) the geocentric 
distance, m & = —27.09 the i?-band magnitude of the Sun (jRussellll916l ). [3 (mag deg -1 ) the 
l inear phase coefficient, and a is the ph ase angle. For 2 IP, we adopt value of Ar = 0.04 
(ILamy et al.ll2005l ; Fernandez et all 1200 ll ) and = 0.035 mag deg -1 which is generally used 
for most comets with unknown 0; however, iFerrfnl (120051 ) argue that a m ean = 0.046 ± 0.013 
may be more appropriate based on the observed phase coefficient for 10 comets. We find 
tn = 1.82 ± 0.05 km, commensu rate with other estimates nucleus size ranging from 1.0 km 
derived by iTancredi et all (120001 ) or ~ 2.1 km value of derived from the snapshot survey of 
Mueller! Jl992f ). assuming A R = 0.04. 



In principle, variations in the photometric brightness derived from the i?-band individual 
images obtained over » 2 — 3 hr intervals during the period of 2005 December 21 through 
22 (Table [2]) can be used to generate a composite light curve yielding estimates of nucleus 
rotational periods. However, our composite light curve of 21P at this epoch, = 2.32 AU 
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exhibits little periodic variation and we are unable to determine a rotational period. 21P 
was also very active during this epoch, with a significant coma contribution which could 
mask subtle variations (±0 .01 mag) in nucleus flux. A rotational period of 9.5 ± 0.2 hr for 
21P has been reported by (ILeibowitz &: Broschlll986l ). 



3.7. Volatiles CO+CO s 



Volati les are frozen as ices or tra pped as gases in amorphous water ice in the nuclei 
of comets (jCaprial 120021 : iPrialnild 120021 ) . Cometary activity occurs when gases are released 
through sublimation or through the exoergic crystallization of amorphous water ice. W ater 
sublimation primarily drives vigorous activity at r h < 5 A U (IMeech fc Svorenl 12004 ) , but 
other volatiles can still be important drivers. At 1.5 AU, iFeaga et al.l (120071 ) found that 
CO2 gas was correl ated with comet 9P /T empers southern hemisphere — the location of the 
strongest dust jet (IFarnham et al.l 120071 ) — yet water vapor was more strongly associated 
with the northern hemisphere. 

Some molecular species (e.g., CO) are released not only from the nucleus (native source), 
but also from within the coma ("distributed" or "extended" source). In the extensively 
studied, bright OC comet C/1995 01 (Hale-Bopp), t here was evidence for dis tributed CO 
sources, possibly arising from icy grains in the coma (IGunnarsson et al.l 120031 ) at large r^, 
whereas at < 1.5 AU the st rong distribute d sources probably arise from the desorption of 
an organic grain c omponent (lCrovisierlll999l). Although icy grains may be the distributed 



source of the CO ( Gunnarsson 



2003 



Bockelee-Morvan et al.l 120021 ). pure- CO ice grains do 



not survive transport into the coma because of their low sublimation temperatures, 20-100 K 
rtPrialnik et "aD l2004h . 



The IRAC 4.5 /zm bandpass encompasses the CO2 i^-band and the CO 1-0 fundamen- 
tal band at 4.26 and 4.67 /zm respectively. The CO and CO2 bands are observed as the 
dominate components at 3-4 fim in ex cess of the dust continu um emission in ISO spectra of 
comets Hale -Bopp and 103P /H artley ( Crovisier et al. 1999al lbl) and in the ambient coma of 
9P/Tempel (IFeaga et al.l 120071 ). After removing the dust continuum, the IRAC images can 
be used to map the CO and CO2 spatial distributions in the comet coma. In principle, the 
emission from the two molecules can be decoupled when we consider that CO and CO2 have 
different lifetimes in the coma and therefore produce different radial distributions. With a 
measured radial distribution, we can also test for evidence of a n extended source of C O or 
CO2. Emission from other prominent volatiles, such as water (IWoodward et al.l 120071 ). lie 
outside the bandpass of the IRAC 4.5 jum filter. 
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The slope of the CO + CO2 emission in the Spitzer images at 4—16 pixels (corresponding 
to projected cometocentric distances of 6.6 x 10 3 km to 2.6 x 10 4 km) can be determined by 
subtracting a scaled IRAC 8.0 /mi image from the IRAC 4.5 /im image (Fig. [2]). The IRAC 
8.0 /im image is dominated by thermal emission from dust, whereas the 4.5 /im image is a 
combination of emission from gas, dust thermal emission, and sunli ght scattered by dust . 



The color of sunlight sc attered by comet dust varies with wavelength (IJewitt fc Meechlll986 



Kolokolova et al.l 120051 ) . and we cannot constrain the spectral slope of the scattered light 
at 4.5 /im with our data. Rather than estimating the 4.5 /im dust flux from our best-fit 
SED ( §3.5p . we subtract a scaled 8.0 /tm image from the 4.5 /tm image to yield an image 
of the gas emission (Fig. [7]). Our best scale factor is 0.025 ± 0.003, which we derived from 
the ratio of the 4.5 /im and 8.0 /im images. We measured the ratio along a 3 pixel wide 
rectangle positioned on the dust tail from 5 to 40 pixels from the nucleus. The error in 
the scale factor is the standard deviation of the pixels in the rectangle. We fit the residual 
surface brightness profile (emission from the gas) in 24° azimuthal steps. The mean surface 
brightness logarithmic slope is —1.00 (median = —1.04, error in the mean = 0.04). The 
derived gas profile is consistent with a long-lived species ejected in a constant outflow from 
the nucleus, i.e, it shows no evidence for an extended source. 

The lifetimes of CO and C0 2 at 2.4 AU are CO = 89 days and C0 2 = 34 days 
( IHuebner et al.lll992l ). For an outflow velocity of 1.0 km s _1 , the characteristic length scales 
are 7.7 x 10 6 km (CO) and 2.9 x 10 6 km (CO2). In our dust subtracted image, there appears 
to be emission out to 60 pixels from the nucleus (~ 1.0 x 10 5 km), but background stars make 
it difficult to measure the logarithmic slope beyond 20 pixels. Both species are long lived on 
scale lengths < 16 pixels even if the outflow velocity is 0.1 km s -1 . From our analysis, both 
CO and CO2 are equally viable candidates for the 4.5 /im excess in the difference image, 
Fig. [3 For comparison, in comet Hale-Bopp (ICrovisier et al.lll999al ). CO2 was the dominant 
emission species at 2.9 AU (comet 21P was observed at 2.4 AU). 

We can measure the intensity of the gaseous emission and derive upper-limits to the 
CO2 and CO contributions. We have attempted to remove the dust from the 4.5 /im im- 
age; however, emission from the nucleus has not been removed. Although we derived an 
effective radius from our optical observations ( §3.61) . the orientation, shape, and tempera- 



ture distribution of the nucleus is unknown at the epoch of the IRAC images. Rather than 
measuring the flux inside a circular aperture centered on the nucleus, we avoid contamina- 
tion from the nucleus by measuring the gas coma inside an annulus with an inner-radius 
of 6 pixels, and an outer-radius of 8 pixels. We convert the measured annular flux den- 
sity, (1.807 ± 0.027) x 10 27 er g s" 1 cm" 2 Hz" 1 , into l i ne em ission fluxes following the IRAC 



prescriptions discussed by the lSpitzer Science Centerl ( 120061 ). where we perform a linear inter- 
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polation of values given in the IRAC spectral response tables^ near the wavelength position of 
our lines, assuming the residual emission within the effective bandpass of the IRAC 4.5 jum 
filter is either entirely from either CO2 (A = 4.26 jum) or CO (A = 4.67 pm). We find 
F C02 = (3.10 ± 0.22) x 10- 14 erg s" 1 cm" 2 , and F co = (2.86 ± 0.20) x 10" 14 erg s" 1 cm 2 . 
The errors include the IRAC absolute calibration uncertainty of 3%, and a dust subtraction 
uncertainty of 11%. 

The gaseous coma has an azimuthally average profile surface brightness slope of p~ l (see 
Table [3]) , thus we can convert the flux at 6-8 pixels to a flux inside a circular aperture of 
projected radius 7 pixels, F tot = 2np 2 F ann /A, where p = 7 pixels, F ann is the flux measured 



in th e annulus, and A is the area of the annulus in square pixels (e.g., iJewitt fc Danielson 



19841 ) . The total emission inside this 7 pixel aperture yields an average column density, 



< N >, in cm 2 , assuming optically thin emission from the expression 



<N>= F, nt 4ttA 



A 



he g band up 2 



(4) 



where A is the observer-comet distance in cm, A(/zm) is the wavelength of the emission, h is 
Planck's constant (erg s), c is the speed of light (/im s _1 ), is the heliocentric distance in AU, 
gband is the solar fluorescence ^-factor at 1 AU in units of s _1 , and p is the projected radius of 
the ap erture, here in units of cm. We have adopted ^-factors given by ICrovisier fc Encrenaz 
rtl983T ) of 2.86 x 10~ 3 s" 1 for the C0 2 v z band, and 2.46 x 10" 4 s" 1 for the CO 1-0 band, 



noting that other cited literature values for these g-factors differ by <10%. From Eqn. HI we 
derive average column densities of: (3.13 ± 0.22) x 10 11 cm -2 , if the emission is solely from 
CO2 and (3.67 ± 0.27) x 10 12 cm -2 , if the emission is solely from CO. 



Bockelee-Morvan et al.l (Il990l ) measured OH radio emission profiles in spectra of several 



comets as a function of r^, including comet 21P, to derive the expansion velocity of that 
species. For most low production rate comets (8.0 x 10 28 < Qoh (molecules s _1 ) < 1.8 x 
10 29 ) the expansion velocities are constant for r^^lA AU. Radio measurements of the 
OH production rates of 21P during the 1985 apparition indicate that Qoh ~ 3.1 x 10 28 
molecules s~ 4 with little change in the der ived expansion velocity, v p = 0.70 ± 0.09 km s _1 
(see Table 4; iBockelee-Morvan et al.lll990l ) over a range of r^. Given the lack of a strong 
rh dependency of v p at distance greater than 1.4 AU, we will adopt 0.70 ± 0.09 km s _1 
for the parent expansion velocity to estimate the expansion velocities of CO2 and CO for 
21P at 2.4 AU from the Sun. With this velocity, we can derive the production rates, Q, in 
molecules s -1 , for each species: 



^http : / / ssc . caltech . edu/irac/spectralresponse . html 
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Q =< N > 2pv x 10 



10 



(5) 



where p is measured in km, and v in km s 1 . The production rates are: Qco 2 — (5.13 ± 
0.75) x 10 25 molecules s _1 , and Q co < (6.01 ±0.89) x 10 26 molecules s _1 . We note that these 
are upper-limits to the production rates for each species as we cannot distinguish between 
the two. 



Mumma et al.l (120001 ) detected CO at 4.67 /im in the coma of comet 21P at = 1.2 AU 
(pre-perihelion) in October 1998 - — the first IR detec tion of this volatile in a Jupiter-family 



comet. However, two w e eks la ter IWeaver et al.l (119991 ) did not detect CO at the same wave- 
lengths. iMumma et al.l (120001 ) suggest that comet 2 IP has a chemically heterogeneous nu- 
cleus, and th at the CO was l imited to a vent that was not illuminated during the observations 
conducted by I Weaver et al.l (119991 ) . This explanation seems reasonable as the comet was near 
perihelion when seasonal variations of insolation are the greatest. Our IRAC images of comet 
21P were obtained at a post-perihelion distance of = 2.4 AU. If the vent responsible for 
the decrease in CO production remains inactive out to 2.4 AU, CO2 would remain the favored 
molecule responsible for the IRAC 4.5 /im excess. 



3.8. Dust Trails 



Determining the cometary dust fraction of the zodiacal dust complex is be coming in- 
creas ingly important. Recent Spitzer discoveries of dust surrounding white dwarfs (IBecklin et al 
20051 ) suggest that mass-loss from cometary bodies (in Kuiper Belt or Oort Cloud a nalogs) or 



tidal disruption of aster oidal objects are feeding their circumstellar dust complexes (jSu et al. 



( jKelley et al 



2008 



200 7 ; iReach et al.l 2005T). Comets lose a large fraction of their mass in dust grains > 100 pm 



Reach. Kellev. & Svkesll2007l : llshiguro et al.ll2007l : ISvkes & Walkerlll992h . 



These grains weakly respo nd to solar radiati on pressure and form distinct dynamical struc- 
tures known as dust trails (jSykes et al.lll988l ). Trails typically appear to follow their parent 
nucleus along the comet's projected orbit, although the "trail" may also lead the comet. 
Dust trails are ubiquitous in Jupiter-family comets, but remain unobserved in OC comets. 
The existence of OC dust trails is inferred from meteor stream studies of Halley-type comets, 
which are derived from the OC. For example, the Halley-type comet 55P/Tempel-Tuttle is 
the parent body of the Leon id meteor stream and at least 80% of JFCs have dust trails 
faeach. Kellev. fc Svkesl 120071 ) . 



The existence of dust trails in comets wi t h per ihelia near 1-2 AU has interesting astrobi- 
ological significance. IReach. Kelley. fc Syked (120071 ) demonstrate that comets have trail mass 
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production rates that are: (1) similar to water mass production rates, and (2) greater than 
the small grain (radii < 10 /im) mass production rates. If the composition of the comet's 
coma approximates the comet interior, then comet nuclei are mostly refractory in nature 
(dust-to-ice mass ratio > 1.0), which limits the amount of water a comet could deliver to 
the surface o f a young te rrestrial planet. Comet 2 IP is the parent body of the Draconid me- 
teor stream (lBeechlll986l). which ca n produce exceptionally active meteor showers (~ 10000 
meteors per hour; I Jenniskend Il995l ) . The outburst activity is strongly associated with years 
when the Earth and 21P closely approach each other, suggesting that the comet is presently 
ejecting meteoroids. 

Thermal radiation from large gr ains in a comet trail is easily detected at IR wave - 
lengths from space-based telescopes (IReach. Kelley. & Sykesl 120071 ; ISykes fc Walker! Il992l ). 
where the emission from warm (~ 100 K) dust peaks. We examined our MIPS 24 /im 
image of 21P (obtained at 2.4 AU) for dust trail emission along the projected velocity vec- 
tor of the comet (Fig. [8]). No trail is detected with a 3a surface brightness upper-limit of 
0.3 MJy sr -1 pixel -1 . Assuming a typical trail grain temperature of « 300 r" 5 = 193 K 
( ISykes fc Walker! Il992l ). the latter surface brightness corresponds to an IR optical depth of 
2.2 x 10~ 9 . O ur upper-limit is larger than the typical comet trail observed in other Spitzer 
MIPS image s (IReach. Kelley. &: Sykesl 120071 ) . thus it is possible that 21P still has a dust trail. 
Miura et al.l (120071 ). using deep i?-band imaging of 2 IP, do not detect the presence of bright- 
ness enhancements at the expected position of the dust trail, and derive a upper-limit to the 
number density of trail grains to be 10~ 10 m -3 . Using similar assumptions (typical grain siz e 
of 1 mm), and assuming a typical trail thickness of ~ 10 4 km (IReach. Kelley. fc Sykesl 120071 ). 
our MIPS optical depth limit sugges ts a number densit y <7x 10~ n m -3 — slightly lower 
than the optically determined value. iMiura et al.l (120071 ) also estimate the dust trail number 
density from meteor shower activities. Their lower-limit (8 x 10~ 15 m -3 ) is about 4 orders 
of magnitude smaller than o ur upper-limit. For c omparison, the trail grain number density 



10 -n m -3 



for 81P/Wild is ~ 10 9 m 3 (llshiguro et al.ll2003l ). and for 67P/Churyumov-Gerasimenko is 



( Kelley et al. 20081 ) . Altogether, the evidence suggests comet 21P has a dust 
trail fainter than the detectable limits of the available observations. 



4. CONCLUSIONS 

We have presented new optical and Spitzer infrared observations of comet 21P/Giacobini- 
Zinner obtained during its 2005 apparition. Analysis of optical imagery indicates that 21P 
was dusty (peak Afp = 131 cm -1 ) and active out to heliocentric distances > 3.3 AU following 
a logarithmic slope with of —2.04. Onset of nucleus activity occurred at a pre-perihelion 
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distance r h ~ 3.80 AU (—375 days pre-perihelion), similar in behavior to that observed in 
the 1991 apparition. The derived average coma colors, V — R = 0.524 ± 0.003,/? — I = 
0.487 ± 0.004 are slightly redder than solar, comparable to colors derived for other Jupiter- 
family comets. Pre-perihelion observations during quiescence yields a nucleus radius of 
1.82 ±0.05 km. 

Spitzer IRAC images obtained at = 2.4 AU, post-perihelion exhibit an extensive 
coma with a prominent dust tail, where excess emission (over the dust continuum) in the 
4.5 /im image at cometocentric distances of ~ 10 4 km likely arises from CO2, although a 
distributed source of CO cannot be discounted. The upper limits to the production rates 
are Q C o 2 < (5.13 ±0.75) x 10 25 molecules s~ l and Q co < (6.01 ±0.89) x 10 26 molecules s _1 . 
The surface brightness of the gas emission is observed to peak along the sun angle, while 
the dust tail peaks near the anti-sunward angle. A search for dust trail emission along the 
projected velocity vector of comet 21P using our MIPS 24 jum image (r h = 2.4 AU), yielded 
no trail (3a surface brightness upper-limit of 0.3 MJy sr _1 pixel -1 ), suggesting that the 
number density of trail particles (typical particle size ~ 1 mm) is < 7 x 10~ n m~ 3 . The 
bolometric albedo of 2 IP derived from the contemporaneous optical and Spitzer observations 
is A{6 = 22°) = 0.11, slightly lower than values derived for other comets at the same phase 
angle. 

This work is based on observations made with the Spitzer Space Telescope, which is 
operated by the Jet Propulsion Laboratory, California Institute of Technology under a con- 
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Table 1. OBSERVATIONAL SUMMARY LOG 



Observation 

Date r A a a FWHM b Filter Integration 

(UT) Telescope (AU) (AU) (cleg) (") Exposures (sec) Sky c 



Pre-Perihelion 



2004 Jun 21 


UH2.2 m 


3.80 


3.21 


13.61 


0.80 


R x 2 


1200 


C 


2004 Jun 22 


UH2.2 m 


3.79 


3.21 


13.76 


0.80 


R x 3 


300 


C 


Post-Perihelion 


2005 Oct 20 


UH2.2 m 


1.76 


1.74 


33.02 


0.60 


R x 11 


30 


P 


2005 Oct 22 


UH2.2 m 


1.78 


1.74 


32.80 


0.76 


V x 2, R x 3 


90 


P 














1x2 






2005 Dec 02 


Spitzer 


2.14 


1.94 


28.19 




MIPS24xl4 


10 




2005 Dec 21 


UH2.2 m 


2.32 


1.72 


22.56 


0.71 


R x 13, I x 2 


200 


P 


2005 Dec 22 


UH2.2 m 


2.32 


1.72 


22.32 


0.65 


V x 3, R x 18 


150 


P 














J x 3 






2005 Dec 31 


Spitzer 


2.40 


1.89 


23.61 




IRAC4.5x60 


30 




2005 Dec 31 


Spitzer 


2.40 


1.89 


23.61 




IRAC8.0x60 


30 




2006 Mar 05 


UH2.2 m 


2.91 


2.24 


16.37 


0.84 


R x 5, I x 1 


500 


c 


2006 Mar 06 


UH2.2 m 


2.92 


2.25 


16.45 


0.73 


V x 1, R x 4 


400 


p 














I x 1 







a Phase angle (Sun-Comet-Earth) or (Sun-Comet-Spitzer). 

b Seeing measured from stellar point sources. 

c Sky shows either cirrus (C) or photometric (P) conditions. 
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Table 2. DERIVED OPTICAL PHOTOMETRIC PARAMETERS 



Observation 



Date 


r h 


i?-band a 


V-R a 


R-I a 


m(l,l,0) 


Afp 


(UT) 


(AU) 


(mag) 


(mag) 


(mag) 


(mag) 


(cm) 




Pre-Perihelion 


2004 Jun 21 


3.80 


21.89±0.04 






17.19 


4.80 


2004 Jun 22 


3.79 


21.95±0.06 






17.24 


4.55 


Post-Perihelion 


2005 Oct 20 


1.76 


15.91±0.01 






12.82 


130.66 


2005 Oct 22 


1.78 


16.05±0.01 


0.53±0.01 


0.49±0.01 


12.89 


124.47 


2005 Dec 21 


2.32 


17.05±0.01 




0.49±0.01 


13.74 


83.00 


2005 Dec 22 


2.32 


17.07±0.01 


0.52±0.01 


0.48±0.01 


13.76 


82.07 


2006 Mar 05 


2.91 


16.37±0.01 




0.48±0.01 


14.80 


39.06 


2006 Mar 06 


2.92 


18.48±0.01 






14.61 


46.50 



a Average color measured in a 3" circular aperture. 
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Table 3. DERIVED MID-INFRARED PHOTOMETRIC PARAMETERS of 21P 



Instrument / 






S.B. 






Wavelength 


rh 


Fx 3, 


Slope, 


k b 




(/im) 


(AU) 


(xlO- 20 W cm" 2 /im- 1 ) 






(xlO™ 20 W cm™ 2 /im" 1 ) 


IRAC 4.5 


2.40 


1.63 ±0.05 


-1.10 ± 


0.07 


0.81 ±0.04 


IRAC 8.0 


2.40 


9.54 ±0.20 


-1.01 ± 


0.06 


4.08 ±0.10 


MIPS 24.0 


2.14 


20.6 ±0.8 


-0.99 ± 


0.03 


6.47 ±0.26 



a F\ measured in a 10" radius (p) aperture. 

b The logarithmic slope of the azimuthally averaged, surface brightness profile. 

C F\ equivalent to a 3741 km radius aperture with the same observing circumstances as 
the 2005 Dec 22 UT optical observations (i.e., corrected for Spitzer-comet distances, dust 
temperature, and activity rate (Afp). 
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Fig. 1. — Selected optical images of comet 21P/Giacobini-Zinner obtained on the UH 2.2 m 
telescope. The fields of view are 180" x 180" in size. North is up and East is to the left in all 
the images, while the stellar point sources appear elongated due to non-sidereal track rates 
of the telescope, (a) 2004 Jun 21, = 3.80 AU pre-perihelion, the comet is marked by the 
ticks in the image; (b) 2005 October 22 = 1.76 AU post-perihelion; (c) 2005 December 21 
rh = 2.32 AU post-perihelion; and (d) 2006 March 05, = 2.91 AU post-perihelion. 
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Fig. 2. — Comet 21P/Giacobini-Zinner Spitzer IRAC images obtained at r h = 2.40 AU 
post-perihelion. The gray-scales range from 0.05 to 0.11 MJy sr _1 for the 4.5 /iin image, and 
1.79 to 3.11 MJy sr" 1 for the 8.0 fxm image. Arrows mark north, east, the projected sun 
angle (©), and the projected velocity of the comet (t>). 



-26- 



3 i m n i i t 1 1 n 1 1 i h 1 1 i i n i 1 1 1 1 1 1 n 1 1 1 1 n i u 1 1 1 1 ■ | i C 



300 '- 24 



200 '- 



_ 100 - 



- 



I 

-100 



-200 - 



— 

-300 " 



n -3 

fa 



v 
E 

O 




h 1 ■ 1 1 1 1 " 1 1 H 



1 



200 100 







^100 
Offest {arcsec) 



-400 -500 



Fig. 3. — Comet 21P/Giacobini-Zinner Spitzer MIPS 24 /xm image obtained at = 2.14 AU 
post-perihelion. The gray-scale ranges from 0.21 to 6.80 MJy sr _1 . Arrows mark north, east, 
the projected sun angle (©), and the projected velocity of the comet (v). 



-27- 



■ 2 



"i — i — i — i — i — i — i — i — i — i — i — i — i — r 



"i — i — i — i — i — r 



14 



— ^ 
■ 

n 



1 6 



J I I I I I I I I I I I I L 



1 



J I I I I L 



Heliocentric Distance [AU] 



Fig. 4. — Comet 21P/Giacobini-Zinner broadband data reduced to unit heliocentric distance, 
r h (AU), geocentric distance, A (AU), and zero phase plotted versus r h . The square symbols 
are for pre-perihelion data and the triangular symbols are for post-perihelion data. The 
dotted horizontal lines represent the likely brightness range for a bare nucleus, calculated 
from our photometry data when no coma was seen (see Table [2]). 
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Fig. 5. — Afp values for comet 21P/Giacobini-Zinner as a function of heliocentric distance, 
rh (AU), for the period between 2004 June through 2006 March. The square symbols are for 
pre-perihelion data and the triangular symbols denote post-perihelion data. 
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Fig. 6. — The spectral energy distribution (SED) derived for comet 21P/Giacobini-Zinner 
from ground-based optical and Spitzer photometry obtained in 2005 December. The smooth 
curve is the best fit curve consisting of a scaled, reddened solar spectrum to represent the 
scattered light at the optical wavelengths, plus a scaled Planck function, to represent the 
thermal emission from the IR wavelengths. 
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Fig. 7. — Comet 21P/Giacobini-Zinner Spitzer IRAC 8.0 /iin image with contours of the 
excess emission (CO/CO2 gas) found in the 4.5 /zm image. The contours are set to 0.05, 0.10, 
and 0.15 MJy sr _1 and the gray-scale range is 1.79-7.08 MJy sr _1 . The surface brightness of 
the dust tail peaks near the anti-sun angle, and the surface brightness of the gas peaks along 
the sun angle. Arrows mark north, east, the projected sun angle (©), and the projected 
velocity of the comet (v). 
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Fig. 8. — A smoothed version of the Spitzer MIPS 24 /im image (contours) and the projected 
orbit of the comet 21P/Giacobini-Zinner (solid black line). The image has been smoothed 
with a Gaussian function, 6.4" FWHM. The contours are logarithmically spaced, with each 
step being a 50% increase in surface brightness, starting at 0.3 MJy sr" 1 , and ending at 
25.9 MJy sr _1 . Boxes mark the areas where we estimated the trail upper-limit from the 
standard deviation of an 11 x 11 pixel aperture. 



